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Q>^ ' Abstract 

\ Recent interest in novel phases in high density QCD motivates the study of high 

■ density supersymmetric QCD (SQCD), where powerful exact results for supersymmetric 
\ gauge theories can be brought to bear in the strongly coupled regime. We begin by 

^ • describing how a chemical potential can be incorporated into a supersymmetric theory 

as a spurion vector superfield. We then study supersymmetric SU{Nc) gauge theories 
with Nf flavors of quarks in the presence of a baryon chemical potential /i, and describe 
pH I the global symmetry breaking patterns at low energy. Our analysis requires < A and 

is thus complementary to the variational approach that has been successful for /i ^ A. 
We find that for Nf < Nc a, modified U{1)b symmetry is preserved, analogous to 

■ the non-supersymmetric 2SC phase, whereas for Nf = Nc there is a critical chemical 
I potential above which the U{1)b is broken, as it is in the non-supersymmetric CFL 

phase. We further analyze the cases with iVc + 1 < Nf < |iVc and find that baryon 
number is broken dynamically for ii > fic- We also give a qualitative description of the 
phases in the 'conformal window', |iVc < Nf < 3Nc, at finite density. 
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1 Introduction 



In the past several years there has been a lot of interest in QCD at high density where novel 
phases of matter are predicted to appear (see the reviews (2 ^ and references therein) . 
High densities imply that the Fermi surface lies at a high energy scale for QCD, which 
means it is weakly coupled. The tools of BCS theory can then be used to demonstrate the 
instability of the Fermi surface and the existence of lower energy condensates which exhibit 
interesting symmetry breaking patterns such as color superconductivity and color-flavor 
locking. These new phases might eventually be observable in compact star dynamics. One 
important consequence of high density QCD is the dynamical breakdown of baryon number 
for three or more quark flavors. 

Because of asymptotic freedom previous treatments of finite density QCD were limited 
to very high density where the physics is perturbative. It would be interesting to find QCD- 
like theories that are calculable even at low densities where the Fermi surface lies within 
the strongly coupled regime. During the last decade, there has been significant progress 
uncovering exact results in supersymmetric gauge theories. In particular, the low energy 
description of supersymmetric-QCD (SQCD) has been described for all numbers of colors 
Nc and fundamental matter fiavors, Nj. (For reviews, see |31I30 

In this paper we use the exact results provided by supersymmetry to study the effects of 
a baryon chemical potential on the vacuum of supersymmetric QCD. The tight restrictions 
imposed by supersymmetry allow us to reduce the set of possible vacua and symmetry 
breaking patterns. A clear limitation of this approach is that the symmetry breaking 
patterns of supersymmetric and non-super symmetric QCD are different even in the absence 
of a chemical potential, so we cannot expect them to be the same in the presence of a 
chemical potential. On the other hand, one can ask a less specific question: is baryon 
number broken dynamically? 

The main result of our study is that indeed this question has the same answer in both 
non-supersymmetric and supersymmetric theories for Nj < |A^c as described below, giving 
additional support to the variational results using the BCS states. We emphasize that 
the dynamical breakdown of U{1)b we have studied differs from the simple Bose-Einstein 
condensation which necessarily occurs in a theory containing scalars when the chemical 
potential exceeds the mass of the scalar. We show explicitly that baryon number is broken 
as a consequence of the strong gauge dynamics at a lower chemical potential. It is also 
intriguing that baryon number can be dynamically broken even if the chemical potential is 
much smaller than the dynamical scale. 

The results for the symmetry breaking patterns for non-supersymmetric QCD in the 
presence of a large chemical potential have been studied in IHI for Nc = 3 and are briefiy 
summarized in Table Note that for Nf < Nc a global U{1) remains unbroken. For 
Nf = 2 this U{1) is a combination of the original baryon number and a diagonal color 
generator. However, for Nf > Nc baryon number is spontaneously broken. Below we will 
compare these results with the possible symmetry breaking patterns in SQCD for Nc > 3 
and various numbers of fiavors. 
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Unbroken Global Symmetry 


1 


U{1)b 




2 


SU{2)l X SU{2)r X U{1)b 


(2SC) 


3 


SU {?>)l+r+c 


(CFL) 


4 


SU{2)v X SU{2)v X SU{2)a 




5 


SU{2)l+r 




6 


SU{3)l+r+c X SU{2)l+r+c X 


U{l)v X C/(1)a 



Table 1: Unbroken subgroup of the original SU {Nf)^ x SU {Nf)R x U{1)b global symmetry 
for QCD with Nc = 3 according to jSHOl- ^ 



We want to study the dynamics of super symmetric QCD in the presence of a chemi- 
cal potential /i which explicitly breaks super symmetry. One main difference between the 
super symmetric and non-supersymmetric gauge theories is of course the presence of scalar 
quarks. If the chemical potential is larger than the scalar mass the squarks immediately 
undergo the standard Bose-Einstein condensation. Since we are interested in the effect of 
strong gauge dynamics, we will add a stabilizing mass to prevent such condensation. Such 
a mass can be either super symmetric or SUSY-breaking. We will study both cases in turn. 

Throughout we will assume that ;U < A so that the chemical potential can be treated 
as a small perturbation compared to the strong supersymmetric gauge dynamics that have 
been well studied. This is a very different regime from the range of validity of the BCS varia- 
tional method in color superconductivity, namely /x » A. In the latter case, the dynamics is 
weakly coupled and the analysis is under control. On the other hand, in our supersymmet- 
ric analysis the dynamics is strongly coupled. The added constraints from supersymmetry 
allow us to draw interesting conclusions about symmetry breaking patterns. Our analysis 
is therefore complementary to the variational method of non-supersymmetric color super- 
conductivity. However, we cannot say anything about gauge non-invariant quantities such 
as (qq). Furthermore, we have assumed that the vacuum is translationally invariant, which 
means we are ignoring the possibility of a crytalline phase |7j. 

Here we briefly summarize our results. We find that baryon number is unbroken when 
Nf < Nc, but that there is a critical chemical potential above which baryon number is 
broken when Nc < Nj < ^Nc- This pattern matches the non-supersymmetric results. 
The dependence of the critical chemical potential on the stabilizing mass can be of two 
qualitatively different types, as illustrated in Figure ^ The situation for Nj = Nc is shown 
in the left diagram, where fic is always less than the stabilizing mass m. The situation for 
more flavors is illustrated in the right diagram, where we see that there is a minimum value 
of the stabilizing mass that allows for a nontrivial critical chemical potential. In all cases 
we require /i < A so that our supersymmetric analysis is valid. 

The outline of the paper is as follows. In the next section we review the formalism for 
constructing an effective Lagrangian in the presence of a chemical potential. In Section 3 

^The result for Nf = 4 in the published version of |3] is incorrect, but has been corrected in the e-print 
version. For Nf = 6 the remaining U{l)v x U{1)a symmetry is a subgroup of SU{6)l x SU{6)r and is thus 
physically distinct from U{1)b- 
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Figure 1: Schematic phase diagrams in the /x-m plane where m is the super symmetric 
mass. The diagram on the left shows the situation for Nf = Nc and the diagram on the 
right shows the situation for Nc + 1 < Nf < 2Nc. The dashed line is the critical chemical 
potential. These diagrams are only valid for fi < A. 

we adapt this formalism to SQCD. In Sections 4-7 we determine the global symmetries of 
SQCD with various numbers of flavors, and in Section 8 we conclude. Appendix A contains 
a simple and explicit example in quantum mechanics, as a check on the method of including 
a chemical potential, while Appendix B reviews the exact results for soft masses in SUSY 
theories. 



2 Relativistic Bose-Einstein Condensation 

Let us begin by reviewing the relativistic formulation of Bose-Einstein condensation for a 
non-supersymmetric scalar field theory. A nice description is given in jBl |[9| . There are two 
purposes to this discussion. One is to show that we can regard the chemical potential as the 
time component of a fictitious gauge field of the U{1)b symmetry at zero temperature. The 
other is to find the criterion for the U{1)b not to be immediately broken in the presence 
of a chemical potential so we can study the dynamical breakdown. 

The partition function in the grand canonical ensemble with nonzero chemical potential 
can be calculated using the following path integral: 

Z = Tre-^(^-^^) 

= C J Vtt^Vtt J V^^V^ J'' '^^^'^''^ 

where M is the time-component of some conserved current. We consider the case of a 
complex scalar field, where the Hamiltonian is 

7^ = vrtyr + V<j)^ ■ V(f> + m'^(l>'<<p (2) 

and the conserved current is = d^(j) — <j)d^<j)^). Thus = i ((/)'!" tt "I" — (/jtt). The integrand 
of the exponent in the path integral can be rewritten 



i7r|i+i7rt^-W(7r,,}i)+M/-(7r,«i) 
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-V0"^ • V<j) - m^^U (3) 

Performing the functional integration over vr and vr^ leads to the following expression for 
the partition function: 

Z = C' J Vcp^Vcf) e~ ■I'o '^■^ / '^'^'^ [(9r+M)0U9r-/i)</'+V<At-V(/.+m20t</.] 

which can be analytically continued to Minkowski space to yield an effective Lagrangian: 



= C' j V(j)^V(j) / /'^'^ (5) 

It is important to note that C^s is not simply C + /iM, since is a function of n in 
addition to (p. Instead, 

£efr = 5>ta^^ + (^ta^0 _ 09^0t) _ (^2 _ ^2)^t^. (6) 

The term linear in /u is the expected fiM contribution. The term quadratic in ^ arises from 
the modification of the conjugate momenta vr = </>^ + ificj). Here we see immediately that 
for > there is a critical chemical potential, \fic\ = m such that for \fi\ > |/ic| the scalar 
potential is unstable and Bose-Einstein condensation takes place. This can also be seen 
explicitly by using the result for the transition temperature for cp'^ theory as in |H|: 

T^ = l{t^'-m'). (7) 

This demonstrates again that condensation occurs for /^^ > m^. 

In Appendix ^ we test this method of adding a chemical potential in quantum me- 
chanics. We do this by comparing the partition function computed directly with the one 
obtained using the effective Lagrangian derived above. The system we consider is a two- 
dimensional harmonic oscillator with an added source term for z-component of angular 
momentum which plays the role of the conserved current. 

The purpose of the preceding exercise was to demonstrate that the inclusion of a chemi- 
cal potential in the Lagrangian amounts to modifying the time derivative dt — > df—ifJ-, which 
is equivalent to the addition of a non-dynamical gauge field that acquires a nonzero vac- 
uum expectation value (VEV) for its time component. Explicitly, the effective Lagrangian 
of Equation © can be rewritten as 

£eff = {Df^cp^D^cP - m^<pU (8) 

where 

D^<P = {d^ - igA^)cl) and {A^) = (^^, 0, 0, . (9) 

It can be shown that the same "covariant derivative" gives the correct effective Lagrangian 
for a theory including fermions and is therefore suitable for use in supersymmetric theories. 
The advantage of this formalism is that the coupling of the chemical potential to the low 
energy degrees of freedom is determined by their U{1)b quantum numbers. 
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3 SQCD with a Baryon Chemical Potential 



We will study supersymmetric QCD which is defined to be an SU{Nc) super Yang- Mills 
theory with Nf flavors of chiral superfield quarks Qi in the fundamental Nc representation 
of the gauge group, along with Nf flavors of chiral "anti-quarks" in the conjugate 
representation. This theory has the anomaly free global symmetry SU{Nf)L x SU{Nf)ji x 
U{1)b X U{l)fi. A frequently appearing quantity is the one-loop beta function coefficient 
bo = 3Nc — Nf. Much is already known about the infrared limit of SQCD, which depends 
on the relative numbers of flavors and colors. For Nf < Nc + 1 the low energy description is 
in terms of composite mesons and baryons, whereas for higher A'^^ it is in terms of mesons 
and dual quarks. 

In this paper we study SQCD with a baryon chemical potential. The regular SQCD 
Lagrangian is defined at a high UV scale, where we also add chemical potential terms 
associated with baryon number for both fermions and bosons. In order to most easily apply 
the known exact results for supersymmetric gauge theories mentioned above 'Wi, we will 
follow the discussion in the previous section and incorporate the added chemical potential 
as a background (fictitious) U{1)b super-gauge field that has received the appropriate 
vacuum expectation value. The resulting UV Lagrangian is 

Cuv = J d^O (Q\e3B^^Qi + Qle-f^^^Qi) + ^ j (W«W" + h.c.) (10) 

where 

{VB) = ea'e{A,) with (^^) = ("^,0,0,0). (11) 

\9B J 

In Equation H10|) the SU (Nc) gauge coupling has been suppressed and the second term 
gives the gauge kinetic terms for the SU{Nc) gauge bosons but not for the U{1)b, which 
remains non-dynamical. 

Note that the VEV in Equation explicitly breaks supersymmetry, so one may worry 
that this invalidates the powerful supersymmetric results, namely that color is confined and 
that the low energy degrees of freedom can be described by mesons, baryons, or dual quark 
superfields. We will therefore consider a chemical potential that is small compared to the 
dynamical scale A. In such a case we can treat the chemical potential as a perturbation to 
the supersymmetric dynamics. 

As we have seen above, the chemical potential /x contributes a tachyonic mass term to 
scalar potentials as in Equation (jEJ. For SQCD with massless quark superfields, this means 
that squark condensation is immediately favored. Because we are interested in studying the 
effect of strong dynamics on condensation and symmetry breaking, we will add a stabilizing 
mass for the squarks that returns the squark VEV to the origin. From this stable UV theory 
we can move on to investigate the symmetry breaking patterns that are triggered in the IR 
description of the theory by strong dynamics. 

We will consider both supersymmetric and SUSY-breaking masses to stabilize the theory 
in the UV. A supersymmetric mass term appears in the superpotential: 

W^mass = miJQ^QJ = TV(mM). (12) 
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In order to preserve the most global symmetry, we take the SUSY mass to have the form 
niij = m6ij, which explicitly breaks the chiral flavor symmetry to the diagonal vector 
subgroup and leaves U{1)b intact. Stability further requires ^ < m. Such supersymmetric 
mass terms are nice because they do not damage the exact results of SQCD. 

If one wishes to maintain the full global symmetry, the stabilizing masses must break 
SUSY. A supersymmetry-breaking soft mass may be added to the Lagrangian: 

Cm^s = -fh^ {Q^Q + Q^Q) , (13) 

where the Q's here represent the scalar components of the respective superfields. Again, 
stability requires ii < fh. Since we rely on many exact results that are only valid when 
the theory is nearly supersymmetric, we will also require m <C A. However, the presence 
of soft masses for the squarks in the UV theory may alter the low energy potential for 
mesons and baryons (or dual quarks). This question was addressed in ^Uj and ^2] using 
spurion arguments (and from a different approach in JHl)- The main observation is that 
scalar soft-masses can be incorporated in a superfield Z whose couplings are determined 
by an anomalous U{1)a symmetry and RG-invariance. Results for softly broken SUSY are 
reviewed in the following sections and in greater detail in Appendix IbI 

We are now ready to investigate the effects of adding stabilizing masses and a chemical 
potential to a UV SQCD theory. In each of the theories below we will first determine 
the symmetry pattern in the absence of a chemical potential, and then check whether 
adding a chemical potential breaks symmetries. We are be particularly interested in com- 
paring U{1)b and its breaking pattern to the non-super symmetric case. Since the IR 
degrees of freedom depend on the number of flavors, we will consider separate cases in 
the following order: Nf = Nc (quantum modified moduli space), Nj < (Affleck-Dine- 
Seiberg superpotential), Nj = N^ + l (confinement without chiral symmetry breaking), and 
Nc + 2 < Nf < 3Nc (dual picture). 

4 Nf = N^ 

When the number of flavors and colors are equal the infrared description is in terms 
of gauge invariant chiral superfields called "mesons" Mij = Q^Qj and "baryons" B = 
^ii---iNjQii' ' ' QiNf ^'^d B = eii-if^^Qii ■ ■ ■ Qij^y Here the suppressed color indices are 

completely antisymmetrized. The mesons transform as {Nf,Nf) under the flavor symme- 
tries, whereas the baryons are flavor singlets. The meson and baryon superfields satisfy a 
constraint that is modified quantum mechanically and can be implemented in a dynamical 
superpotential with a Lagrange multiplier superfield X: 

W = jm^ (detM -BB- A2^=) , (14) 

where A is the scale of the SU{Nc) theory. The equation of motion for X enforces the 
constraint. 

We will now add stabilizing masses. In the following subsection we will add supersym- 
metric masses. In Section [4.21 we will add soft masses and discover that we need to choose 
a vacuum to expand around. In 14.2. II and 14.2.'^ we will consider two such vacua in turn. 
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4.1 Supersymmetric Mass 

As discussed above, the theory with massless squarks in the UV is immediately unstable 
upon the introduction of a chemical potential. We will first consider supersymmetric mass 
to stabilize the theory in the UV. When the supersymmetric mass term in Equation (|12|) 
is added to the superpotential (fTl)) the F-term equations force the baryons to vanish, but 
yield the VEVs 

(M,,) = {m-% {a'o detm)'^""^ and (X^^) = (15) 

The VEV for X then manifests itself as a mass term for the baryons. The low energy 
Lagrangian for the baryons is thus 

C = J d^ed^9j^iB^e''^^^B + #e-^^^^^) + J d^j^BB (16) 

where cq is an unknown coefficient in the Kahler potential which is estimated in naive 
dimensional analysis to be 0(1). We've shown that at zero chemical potential U{1)b is 
conserved and the chiral symmetry is broken to its vector subgroup. 

We will now show that this chang chemical potential is added. After canonically 

normalizing B and B and inserting the VEV of Vb from Equation Hll|) to incorporate the 
chemical potential, the squared mass of the baryons is 

ml^^ = mVco - N^,^?, (17) 

which means there is a critical chemical potential 

2 

such that for n > Hc the baryon mass-squared becomes negative, drawing the baryon VEVs 
away from the origin. We cannot determine exactly where the baryons stabilize, but due 
to the presence of the SUSY masses we know that the UV theory is stable, which means 
the baryon VEVs should be of order A. Thus we conclude that as long as cqN^ > 1 there 
are values of such that /Xc < /x < m where our approximations are valid and the U{1)b 
symmetry is spontaneously broken. This is depicted in the left diagram of Figure ^ 

4.2 SUSY-Breaking Mass 

Instead of adding supersymmetric mass terms, we now consider the addition of soft SUSY- 
breaking mass terms to stabilize the UV theory. In order to determine the potential for 
the low energy degrees of freedom due to SUSY breaking, we will adopt a dimensionless 
convention for the fields jEj 

(19) 

H H 
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The constraint then takes the convenient form 



det M - BB = I. (20) 

The fields M, B and B are dimensionless and are uncharged under both the anomalous 
U{1)a symmetry described in Appendix^ and the anomaly free U{1)r. This means we 
cannot use these symmetries to constrain the theory and cannot determine the IR soft 
masses exactly. However, we will still be able to maintain some control over the sign of the 
baryon masses and thus rule out some symmetry preserving vacua due to their instability. 

The flavor and U{1)b symmetries constrain the Kahler potential to be a some function 
K of real and invariant combinations of the fields 



CD J (fe(fe I X K (^B^e'^'^^'^B, i^e'^'^^^i, (-B^ + h.c), trM"^M, . . 



(211 



where we've only written lowest order terms in the fields. Here I is a U{1)a and RG- 
invariant superfield defined in Equation (|82|) . The point of Equation (|21|) is that /, which 
contains the soft mass fin? as its 9^9^ component, only appears as an overall multiplicative 
factor, while the background U{1)b gauge field which contains the chemical potential /i 

couples to the baryons as expected. Note that the BB term cannot be forbidden by the 
symmetries, and prevents us from completely determining the baryon soft masses. In later 
cases such terms will be disallowed by U{1)a invariance. 

In what follows we will assume that the SUSY-breaking effects are sufficiently small 
compared to the dynamical scale A such that the quantum modified constraint, Equa- 
tion (|2()j) . is strictly enforced. We will therefore restrict ourselves to the super symmetric 
fiat directions while enforcing the constraint by hand. We then minimize the resulting 
potential including the soft-breaking terms. 

Since the origin of moduli space is excluded by the constraint, some of the global 
symmetries are necessarily broken by the strong dynamics even in the absence of a chemical 
potential. We will therefore expand the low energy Kahler potential around points that 
possess the highest degree of remaining symmetry and see whether they remain stable in the 
presence of a chemical potential. In the following subsections we will analyze the two points 
of "maximal" symmetry. We cannot determine which point gives the true vacuum of softly- 
broken super symmetric gauge theory. However, the first point resembles the dynamics of 
non-supersymmetric theory the most and hence is a good candidate for the true vacuum. 
In that case we can show that a large enough chemical potential destabilizes the vacuum, 
leading to the spontaneous breakdown of U{1)b- 



4.2.1 Vacuum with Chiral Symmetry Breaking 

In terms of the hatted fields, a general vacuum that conserves baryon number has {B) = 
{B) = and (det M) = 1. Among them, the solution to the constraint (Mij) = 6ij preserves 
the largest symmetry. At this vacuum of the supersymmetric theory, the SU{Nf) x SU{Nf) 
chiral symmetry is broken to the diagonal SU{Nf), while U{1)b is unbroken. Therefore 
this point resembles strongly the dynamics of the non-supersymmetric theory, and hence it 
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is reasonable to expect that the soft supersymmetry breaking effects make this point stable 
by giving the fluctuations positive squared masses. Indeed, we will show below that this 
is not unreasonable. Furthermore, since this pattern is similar to that achieved by adding 
super symmetric masses it is interesting to compare the results for the two types of masses 
once a chemical potential is added. 

We expand the meson field as M = 5ij + ^ij. Following ^j, we can now use Equa- 
tion (|2()j) to eliminate one of the degrees of freedom. It is convenient to express the trace 
of $ in terms of the other fields, 

tr$ = BB + ^tr(<I>'^) + higher terms in tr<I>, (22) 

where = — Sij/Nf is the traceless component of 

The Kahler potential in Equation ()21() can now be expressed in terms of the small 

excitations B and B and expanded in a power series. For example, substituting into the 
meson kinetic term in the Kahler potential we get 

trM^iQ" = Nf + tr($''f$') + ^ (tr($'2) + h.c.) + (^BB + h.c.) + . . . 

1 2 

= Nf + -tr + $'t) + (bB + h.c.) + . . . (23) 

The soft masses will ultimately be determined by inserting the SUSY breaking VEVs 
for Vb and / (Equations and (|83l respectively') and performing the integral in the 
Kahler potential. The mesons will receive a soft mass only from the SUSY breaking in /. 
From Equation H23|) we see that at lowest order in the fields only the real part of receives 
a mass while the imaginary part remains massless. However, since we are expanding the 
meson field around a VEV of order the dynamical scale, masses that come from inserting 
the VEV into the higher order terms in Equation (|21|) will not be suppressed. 

We can use Goldstone's Theorem to show that the imaginary part of remains massless 
to all orders in the fields. The Kahler potential of Equation (j21j) can be expanded in powers 
of the excitations'^ 



£ D ^ d^ed^e I X 

[c2BB + h.c.) + catr^)'^^' + C4 (tr^'^ + h.c.) ] . (24) 



The coefficients Cj are derivatives of the function K and include contributions from operators 
of all orders in the meson field. The coefficients ci and C3 should be absorbed into the 
definition of the fields in order to have canonically normalized kinetic terms. Performing 
the superspace integration, the mass matrix for each <I>^j in the basis ($', <&'^) is 

ml,=m'( ] ''{'A. (25) 

V C4/C3 1 



*We are assuming that the strong dynamics respect a Z2 symmetry that interchanges B and B. 
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Since the VEV of M breaks the chiral flavor symmetry SU{N)l x SU{N)ji to the diagonal 
SU{N)l^^_ there must be Nj — 1 massless Goldsone Bosons. Since there are A^^^ — 1 <!>' 
fields this can only work out if all the matrices in Equation ()25() have a zero eigenvalue. 
This sets the ratio C4/C3 = 1 (though -1 will do as well) and the masses are 

"^Re*' = "^L*' = (26) 

Therefore, the fluctuations in the meson degrees freedom indeed have positive mass-squared 
for their real parts and zero mass for the Nambu-Goldstone bosons, and hence the vacuum 
is stable against these fluctuations. Since the mesons are uncharged under baryon number 
this is unchanged in the presence of a chemical potential. 

On the other hand, the situation with the baryonic degrees freedom is less clear. In J2j 
the meson masses were determined as above and it was noted that since we do not have 
control over the coefficients of the BB-type mass, the baryon masses cannot be explicitly 
determined. Nonetheless, we can still learn about the stability of the potential. After 
canonically normalizing and performing the 9'^ integral, the baryon mass matrix in the 

basis {B B ) is 

""^ " [ 'fm^ m2 A ^2^2 j (27) 

The diagonal mass receives a positive contribution from the stabilizing soft mass fh? and 
a negative contribution from the chemical potential which is enhanced by a factor of A'^. 
Recall that the requirement for a stable theory is > fi^. 

In order to determine whether a chemical potential induces a phase transition we first 
must check whether baryon number is conserved at // = 0. In the absence of the chemical 
potential, the eigenvalues of the mass matrix Eq. (|27j) are (1 it ^ )m^. Therefore the 

vacuum B = B = {)\s stable if the off-diagonal element is smaller than the diagonal element, 
|c2| < |ci|. We cannot rigorously justify that this is the case, though it is a reasonable 
assumption.^ Henceforth we assume that this is the case. Then the vacuum is stable 
against baryonic fluctuations and the dynamics is precisely that of the non-super symmetric 
theory: dynamical chiral symmetry breaking with massless Nambu-Goldstone bosons and 
massive baryons. The main difference, however, is that the baryons are not as heavy as A 
but rather as light as m <C A (within the validity of our approximations). 

Now we consider a finite chemical potential. As seen from the mass matrix Eq. H27() . 
there is a critical value of the chemical potential that makes one of the eigenvalues negative: 



/^c = ^ ( 1 



m^. (28) 



Without loss of generality, we can always perform a U{\)b rotation to make C2/C1 < 0. 
Then the direction of the instability \s B = B. The fields roll down the potential and 

■''For Nc — 2 the SU{2) x SU{2) flavor symmetry is enlarged to 5*0(4), which prohibits the BB terms. 
However, the SO (4) is broken by the chemical potential which would be expected to generate a BB term 
with a coefficient of order fi/A, which is less than the 0(1) coefficient of the B^ B terms. Thus \c2\ < \ci \ is 
the natural expectation for = 2. 
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dynamically break the U{1)b symmetry. On the other hand, \B\, \B\ ^ 1 corresponds to 
the semi-classical regime where 





( " 










Q = Q = 




J- 


M = 














[ 





B = B = v^^. (29) 



We know the potential is stable in this regime because the dynamics is correctly described 
in terms of the quark degrees of freedom and fh > ^ assures the stability. Therefore, 
the baryon fields should stabilize somewhere around B = B ^ 1. The precise vacuum 
expectation values depend on the exact form of the Kahler potential and hence cannot be 
worked out. However the existence of a stable vacuum is obvious from this discussion. 

In summary, up to the assumption of |c2| < |ci|, the supersymmetric vacuum with chiral 
symmetry breaking and unbroken U (1) b is stable in the presence of the soft supersymmetry 
breaking.^ There is a critical chemical potential of the order of supersymmetry breaking 
beyond which U{1)b breaks dynamically. The order parameters are expected to all be on 
the order of the dynamical scale and hence the theory is strongly coupled. Nonetheless the 
breakdown of U{1)b and the presence of a stable vacuum is guaranteed. 

This same conclusion holds even if we were to expand around more general baryon num- 
ber conserving vacua, since the crucial piece of information we need is only the instability 
of the baryon directions, which follows directly from the diagonal soft mass and chemical 
potential contributions to the baryon mass terms. Since this is independent of the precise 
meson VEV, we conclude that (modulo one assumption) all baryon conserving vacua are 
unstable and thus that baryon number is not a global symmetry in the presence of a chem- 
ical potential that is larger than fic- As in the case with the supersymmetric mass, the 
critical chemical potential is again of order m/Nc- The breaking pattern also agrees with 
the result for non-super symmetric QCD. 



4.2.2 Vacuum with Baryon Number Breaking 

In the supersymmetric limit, there is another vacuum with a large unbroken symmetry, 
{M) = and (BB) = -1. It preserves the fuh SU{Nf) x SU{Nf) chiral symmetry but 
breaks baryon number symmetry. We already see that this point is of lesser interest to us 
since we are focusing on the C/(l)B-breaking signature induced by the chemical potential. 

Again we would like to determine whether this point is stable in the presence of soft 
supersymmetry breaking. The dynamical degrees of freedom around this point are the full 
meson matrix M and another chiral superfield b: B = e^, B = e~*. Since b ^ b + i6 under 
a U{1)b transformation, the Kahler potential depends only on 6 -|- 6*. In the absence of a 
chemical potential, the Kahler potential for mesons takes the form 

£d J d'^eif{b + b*)TiM^M. (30) 

''If \c2\ > |ci|, baryon number is already broken in the absence of a chemical potential. In this case a 
chemical potential does not change the dynamics qualitatively but it does cause the L''(l)s-breaking to be 
"stronger". This will cause the critical temperature at which U{1)b is restored to grow with /i as one would 
naively expect (see Equation 0. 
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The positivity of the kinetic term for meson fields requires /(O) > 0, while the [In/]g2g2 = 
—2^w? gives the mesons a positive mass squared. Therefore, the meson directions are 
stable around this point. 

Concerning the direction 6, the leading term in the Kahler potential is 

j <feh{h + h*f, (31) 

which gives mass squared fh? to the real part while keeping the imaginary part a massless 
Nambu-Goldstone boson. A linear term J d'^9I{b + b*) is forbidden by charge conjugation 
invariance Q ^ Q that requires the invariance of the effective Lagrangian under b —b. 
Therefore this point is stable even in the presence of the soft supersymmetry breaking. 

Because baryon number is already broken in the absence of the chemical potential, the 
finite chemical potential can change the physics qualitatively only by breaking the chiral 
symmetry. However, it is the breaking of baryon number that prevents us from calculating 
whether this is indeed the case. One possibility is that the meson direction is destabilized 
by an operator of the type 

j d9^d9^ cl (^B^e^^^'^Bj (^tiM^ D cN^iJ,'^\Bo\'^tTM^ M (32) 

which could give an additional contribution of either sign to the meson mass squared. Thus 
we cannot determine the fate of the chiral symmetry. In any case, this vacuum is of limited 
interest to us. 



5 Nf<Nc 

When Nf < N^. the non-supcrsymmetric results leave an unbroken global U{1). As we will 
show below this is true in SQCD as well. In this case the IR degrees of freedom consist 
only of the meson superfield M = QQ since there are not enough quark flavors to form 
color-singlet baryons. The SQCD dynamics generates the so-called AfHeck-Dine-Seiberg 
superpotential, 

\3+2Nf/{Nc~Nf) 

W = (Nc- Nf)- , . (33) 

^ ■^^(dctM)V(^c-N/) ^ ^ 

Because of this potential, the meson field is driven away from the origin, leading to a 
run-away theory with no vacuum. 

However, this runaway behavior is stopped by the masses we add in the UV. Let us 
begin with a super symmetric mass and assume that this mass is large, m > A. Adding this 
mass to the superpotential and solving for the meson VEV we get the same formula as in 
the previous section 

/ ^ N l-Nf/N^ / A \ ^/^'^ o 

{Mij) = {m-%(A'°detm) " ~ (^-j A^. (34) 

Note that the VEV is smaller than A and therefore the meson is indeed the appropriate 
degree of freedom to describe the IR vacuum. We can immediately see that baryon number 
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is unbroken because the IR degrees of freedom are not charged under U{1)b- For this case 
adding a chemical potential has no effect. 

We will now consider adding small masses m < A. In this case the supersymmetric 
masses and the soft masses yield similar results so we will treat them together. When the 
stabilizing mass is small the mesons runs far away from the origin before it is stabilized at 
VEVs much larger than A. This region of moduli space is described in terms of quarks. 
The large meson amplitude corresponds to the fiat direction 



/ Vl 



Q = Q 



\ 



VNf 



(35) 



V ••• / 

Rewriting the ADS potential in terms of quarks and adding soft or supersymmetric 

masses and the tachyonic mass due to the chemical potential we can minimize and solve 
for Vl, . . . ,VNp- For example, when soft masses are added such that < rfi^ — (J^, we get 



(Nc-Nf)/2Nc 



Vl 



VNf 



A 



— 



(36) 



The form of the VEV is slightly different when a supersymmetric mass is used, but vi = 
. . . = vjsj^ is still satisfied. Therefore this theory dynamically breaks the chiral symmetry 
SU{Nf) X SU{Nf) SU{Nf). However, a U{l)j^ remains unbroken. The new baryon 
symmetry is the simultaneous transformation under the original U{1)b and a U{iy , where 
U{iy is a subgroup of the gauge symmetry embedded as SU{Nc) D SU{Nf ) x SU{Nc — 
Nf) X U{iy . This situation is analogous to the non-supersymmetric case of Nf = 2 and 
A^c = 3 (the so-called 2SC phase) where the unbroken U(l) is a combination of baryon 
number and a broken color generator. 



6 Nf = Nc^l 

When Nf = A^c + 1 the IR theory is again described in terms of mesons Mij = QfQ'j 

and baryons Bi = €ij^...j^^Qj^ ■■■QjNc = ^m-jNcQji '"QjNrJ symmetry 

transformation properties are different. Choosing the quarks to have baryon number 1, the 
fields M, B and B transform as {Nf,Nf)o, {Nf,l)i\f^ and {l,Nf)-N^ respectively, under 
the global symmetry SU{Nf) x SU{Nf) x U{1)b- 

We summarize our results before presenting the details. Recall that for A^^ = 4 in 
non-supersymmetric QCD baryon number is broken in the presence of a high chemical 
potential. For SQCD, adding a supersymmetric mass in a specified range will again lead 
to a critical chemical potential above which baryon number is broken, in agreement with 
the QCD result. Adding UV soft masses yields only partial agreement. If we choose to 
stay along a flat direction that breaks the chiral symmetry to a vector symmetry we indeed 
find that the existence of a critical chemical potential is plausible. However, there is also 
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another flat direction along which the chiral symmetry is unbroken which has tachyonic 
masses for the baryons, thus breaking U{1)b even for = 0. 



6.1 Supersymmetric Mass 

As in the previous case we wih stabihze the UV theory first with supersymmetric masses 
m and then with soft SUSY-breaking masses rh. 

The low energy superpotential for baryons and mesons, including the supersymmetric 
quark masses, is 

W = -^{detM - BiMijBj) + Tr(mM). (37) 
The baryon VEVs vanish but the mesons acquire a nonzero VEV: 

{Mij) = {m-% (a^o det rnj . (38) 

The meson VEV yields a mass for the baryons / cP9{{M) / A^^')BB. Including the chem- 
ical potential contribution, the mass-squared for the canonically normalized baryon and 
antibaryon flelds is 

<^ = 1(^) A^-iV,V (39) 

where again cq is an unknown 0{1) coefficient in the Kahler potential, and we have taken 
the supersymmetric mass rriij = m6ij in order to preserve as much global symmetry as 
possible. 

Here again there is a critical chemical potential above which the baryons become 
unstable and break U{1)b- 

, 1 fm\y"-. 



However, this formula will only be valid as long as /j,'^ < A^ (near-SUSY limit) and //^ < 
(UV stability). These two constraints yield the requirements 

(coiV,2)^^/(2-2^^) < ^ < (coA^2)JV./2^ (41) 

These can be satisfied as long as cqN^ > 1, since (for Nc > 1) the lower bound is a decreasing 
function of cqN^ and the upper bound is an increasing function. This is the same restriction 
on Co that we found in the case where Nf = Nc- We again conclude that it is very likely 
that baryon number is broken by a sufficiently large chemical potential. There are, however, 
some diff'erences compared to the Nf = Nc result. The critical chemical potential is not 
0{m) but is a combination of m and A. Also note that the UV mass is bounded from below 
in order for our result to hold, as shown in the right diagram of Figure ^ 

6.2 SUSY-Breaking Mass 

Now we consider a different situation where m is set equal to zero in Equation ()37() and 
instead we add universal soft SUSY-breaking masses fh for the squarks in the UV theory. 
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Using the results of ^Oj reviewed in Appendix^ the soft-masses for the squarks result in 
the following soft masses for the mesons and baryons: 



'mil 



2Nr 



2Nr - 1 



-m 



and 



N. 



2Nr - 1 



c ~2 

m 



(42) 



Notice that for Nc > 3 the meson mass-squared is positive, but the baryons are tachyonic. 
Here we have also included the contribution from the chemical potential, which simply 
follows from baryon number. 

When we canonically normalize the baryon and meson superfields the unknown con- 
stants Cj^j ^ ^ reappear as relative coefficients between the terms of the superpotential: 

detM' ; 



W = ci 



(43) 



where M' = M/A^ and B' = B/A^"^^ are dimension one fields, and ci = Cj^,/'-"" and 
C2 = (cmcbc^)^^^"^ ■ This, along with the soft masses of Equation H42() . gives rise to the 
following scalar potential 

2 



Nf/2 



ViR 



m, 



ci 



C2B[B,^ 



+ \c2Mi^B'A +\c2B[M'ij 



+ E 



M'- 



m 



B 



B 



(44) 



^3 



where m'^j is the ijth cofactor of M'^j. In the supersymmetric limit this potential has two 
flat direction with different symmetry breaking patterns. We will now analyze these flat 
directions separately. 

The first flat direction has VEVs of the form: 



(45) 



Along this flat direction the relationship between the VEVs is constrained to be cia'^'' = 
C2A^''^'^b^ by the first term in Equation (|44() . The introduction of the soft masses has the 
potential to destabilize this fiat direction. The potential along the flat direction for \a\ < A 
is given by 
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ViR = N, 



2Nr 



2Nr-l 



~ 2| |2 

m a 



N, 



2Nr 



ci a 



C2 



AN, 



(46) 



On the other hand, at large amplitudes the quark description is valid. The flat direction 
above corresponds to the quark VEVs 



/ 



Q = Q 



\ 




(47) 
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and hence, for \a\ > A, 



ViR 



fi )\a\ 



(48) 



Clearly Vjr = at the origin. The potential goes up as |ap first, but then goes down as 
— |a|''^'' at a larger amplitude. We will have a vacuum with energy lower than the origin if 
the potential in Equation (|46j) becomes negative within its range of applicability, namely 
\a\ < A. Beyond this region the potential goes up linearly with \a\ as in Eq. ()48() . 
When /U = 0, Vjr will be negative when 



l/(Nc-2) 



< A <^ 



(49) 



where the upper-bound is the requirement for Equation (|46|1 to be valid. This condition 
is somewhat restrictive, though by no means impossible to satisfy. For A'^c = 3 it requires 
|ci| > A^c|c2| with \a\ near its upper bound. 

Nonzero /U helps make Vin more negative. In particular, there is a critical chemical 
potential given by 



N, 



1) 



above which Vjr goes negative, as long as 

A^c(A^c - : 



A^. 



Nc-2 



m 



(50) 
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C2 
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N^{2N, -l) + {N,-2) 



l/(iVc-2) 



< 



\a\ 

T 



< 1. 



(51) 



The lower bound comes from requiring fj, < m, whereas the upper bound is again required 
for Equation (|4Hj) to be valid. The extra factors of Nc on the left-hand side make the 
requirements on ci and C2 less stringent. However, as Nc grows the allowed range for \a\ 
decreases. If indeed such a point is the global minimum, baryon number will be broken, and 
the chiral flavor symmetries will be broken to the diagonal SU{Nf — l)v vector symmetry 
once n exceeds fic- 

We will now analyze the second flat direction. When M = a SUSY vacuum requires at 
least one of Bi or Bi to also vanish. For example, take Bi = for all i. The positive meson 
mass-squared will tend to keep the meson VEV zero, but the negative baryon mass-squared 
will draw the B VEVs away from the origin. Specifically, flavor symmetry allows us to 
write the VEVs as 

f b\ ( b\ 



and B 



B 







(52) 



which will be pulled in the direction |6p ^ 0; I^P 
flavor symmetry is broken to SU{N(.) x SU{Nf). 



= or vice versa. In the former case the 
Once the amplitude is larger than A, this 
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direction corresponds to 



Q 



( 



V 





J 



Q = 0. 



(53) 



^2^|^|2/Afc large amplitudes. Thus the 



Again the potential goes up as Vir = Nc{m 
vacuum of the theory will settle at an intermediate scale where we lose calculability. 

Since we did not truly minimize the potential, the two symmetry breaking patterns 
above should be taken as candidates. One may hope that once the theory begins to roll 
down one of these directions it will 'try to maintain as much of its symmetry' and stay 
along that direction. However, it is possible that the true minimum lies elsewhere. 

In summary, when Nf = Nc + 1 and soft masses are added, the origin of moduli space 
is unstable and thus some of the SQCD global symmetries must be broken. We found 
two candidate directions. In the first case the theory is not unstable right at the origin 
but becomes so near the origin as long as Nc is not too big. In this direction the chiral 
symmetries are broken to the vector subgroup, and the meson VEVs as in Equation (|45|) 
show a color-flavor locked pattern for Nc of the flavors. However, in the second direction, 
which is unstable for any number of colors, some of the chiral symmetries remain unbroken. 
It appears that there is a first order phase transition between the two directions as fj, is 
increased. 

In both cases baryon number is broken dynamically, consistent with the result in non- 
supersymmetric QCD. Both of the candidate symmetry breaking patterns differ from the 
previous results for non-supersymmetric QCD as shown in Tabled 



7 Nc-\-2< Nf <3Nc- Dual Region 

In this section we would like to analyze SQCD when Nc + 2 < Nf < 3Nc- Seiberg has 
shown that the original SU{Nc) "electric" theory is dual to a different "magnetic" gauge 
theory. The dual theory consists of an SU{Nc) gauge theory, where Nc = Nf — Nc, with 
Nf flavors of magnetic quarks qi and qi along with a gauge singlet Mij and superpotential 
W = ^qiMijqjJ The two theories are obviously different in the UV, but they are dual in 
the sense that they describe the same IR physics. 

When Nc + 2 < Nf < ^Nc, the so-called free magnetic phase, the dual description is 
particularly useful because it is IR free. In that case, the theory of dual quarks and the 
meson is calculable in the IR and our goal is clear: we want to determine whether the mass 
squared of the dual quarks is negative when a chemical potential is added. 

When ^Nc < Nf < 3Nc both the electric and the magnetic theories are asymptotically 
free but flow to an IR fixed point. Our analysis hinges on the ability to identify the correct 
low energy degrees of freedom, their masses, and their charge under U{1)b- Once these are 
specified the effect of the chemical potential is determined by gauge invariance and we can 

^Here we use k instead of the conventional fi in order to distinguish it from the chemical potential. 
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determine the stability of the origin. However, when the IR physics is strongly coupled the 
quarks (electric or magnetic) may no longer be the correct degrees of freedom. Therefore 
the fixed point value of the gauge coupling is crucial in determining the reliability of our 
result. As Nj is increased within this range the fixed point coupling for the electric theory 
decreases whereas the fixed point coupling of the dual magnetic theory increases. Once 
Nf > 3Nc the electric theory becomes IR-free and the magnetic theory is asymptotically 
free. Therefore, when |A^c < < 2A?c the magnetic theory is more reliable, whereas 
for 2Nc < Nf < 3Nc it is the electric theory wc should be using. Near Nf = 2Nc both 
theories are strongly coupled, so our analysis breaks down in that region. Unfortunately we 
cannot determine the fixed point couplings beyond the qualitative discussion above. Since 
all we can do is calculate the mass squared of either the electric or magnetic quarks we will 
therefore trust the result only in the two edges of the conformal window where one of the 
descriptions is known to be weakly coupled. 

Before we discuss our results we will further comment on calculating the critical chemical 
potential in the conformal window. Below we will use exact supersymmetric results to 
calculate the masses of the IR degrees of freedom, cither electric or magnetic. The mass 
calculated in this manner is the effective mass at the scale A, where the theory becomes 
conformal. If one would like to calculate the masses at a lower scale uj, the aproach to the 
fixed point must be considered. The wavefunction renormalization of the quarks is given 

where 7* = bo/Nf (or bo/Nf in the magnetic case) is determined by the anomaly free R- 
charge. This decrease in Z leads to an enhancement in the effective IR masses ruphys = rajZ. 
One may worry that the increase in the stabilizing mass might make the breaking of U (1)^ 
by a chemical potential more difficult. However, this is not the case. We may write the 
effective potential for the squrk at the scale lo ^ q t^s, 

V,ff{q)=m\q)\q\^. (55) 

If the effective mass m'^{q) is negative for some scale (not necesarily the IR) the potential 
becomes negative there and a non-trivial minimum exists. Since the supersymmetric mass 
is only enhanced below A, our best chance of achieving an instability is for a scale of order A. 

The STimmary of our result for this case is as follows. When we add supersymmetric 
masses and focus on the free magnetic phase we find that a large enough chemical potential 
will make the dual quark masses negative and break U{1)b, provided a constraint on k/A 
is satisfied. Notice that this phase does not exist for A^c = 3 so a comparison with QCD 
results cannot be made. 

For Nf = ^Nc there is no critical chemical potential within our approximations. For 
^Nc < Nf < 2 Ac a critical chemical potential may exist if a certain restriction on k/A is 
satisfied. However the naive expectation of A; ~ A does not satisfy the restriction. This 
is to be compared with the QCD result that U{1)b is broken for Nf = 5 which is at the 
bottom of the conformal window (Nf = 1.66Nc)- 

If we extrapolate our results to Nf = 2Nc we find that the critical chemical potential 
needed to destabilize the magnetic quarks approaches the electric quark mass. This shows 
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that the theory maintains its self duaUty for Nf = 2Nc when a chemical potential is added, 
and might indicate that U{1)b can no longer be broken dynamically in either description. 
However, this is the region where we have the least control since both descriptions of the 
theory are strongly coupled. 

Extrapolating further to 2Nc < Nf < 3Nc we again find that the mass squared for 
the magnetic squarks may become negative with a sufficient chemical potential (without 
requirements on k), however in this region the electric quarks are more reliable. By our 
construction the electric quarks do not condense dynamically and therefore U{1)b is not 
broken for this case either. 

When soft-masses are added in the UV theory we have control over the IR masses 
only in the free magnetic phase. There baryon-number is broken even without a chemical 
potential. Even though a finite chemical potential leads to a "stronger" breaking of baryon 
number, it does not lead to a phase transition like that found in other cases. 

In the conformal window, in addition to the difficulties described above, the IR soft 
masses are uncalculable. Qualitatively we expect that near Nf = ^Nc baryon number is 
dynamically broken by a chemical potential and near A'^^ = 3Nc it is not. This will be 
argued in Section 17.2.21 



7.1 Supersymmetric Mass 

First we will consider adding a supersymmetric mass for all the quarks in the original 
theory in the form of a superpotential VFmass = rriijQiQj, with rriij = mSij as before, which 
preserves the vector flavor symmetry. In the dual theory, we thus have a superpotential 

Wdual = p^M,jqj + Tr (mM). (56) 

The equations of motion for M yield qiijj = rriji, which does not admit any solutions 
because qiqj is at most a rank N^ matrix, whereas niij is chosen to have rank Nf > A^^. 
However, since we know that (M) 7^ in the electric theory, it suggests that the same must 
be true in the magnetic theory. Then (M) becomes a mass term for the dual quarks q, so 
we can integrate them out, leaving a pure SU{Nc) Yang-Mills theory with an additional 
gauge singlet M. The Yang-Mills theory undergoes gaugino condensation, which gives a 

superpotential term Wxx = NcA^j^. The usual matching conditions give Af^ = A''" det(M), 
so together with the superpotential for M we find an effective superpotential 

Wes = Nc I A^^^-^f det(M) I + Tr(mM). (57) 

Minimizing with respect to Mij gives an expression for rriij which, together with the duality 

relation A^oU'o = {-l)^''k'^f, can be inverted to give (Mij) = {m~^)ij ( A^o det m) , 
which is consistent with the assumption that (M) 7^ and matches the result in the 
electric theory. 
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7.1.1 Free Magnetic Phase 

In the free magnetic phase we can estimate the size of k as follows. In the duality relation, 
A is the scale where we expect the electric theory to become strongly coupled, whereas A is 
the high scale where the IR-free magnetic theory is strong. Thus it is reasonable to identify 
these two scales, so that the theory is described by the UV-free electric theory that flows 
down and becomes strong around A = A where the IR-free magnetic theory takes over at 
strong coupling and subsequently flows to weak coupling in the IR. With this identification, 
k is also naturally the same size as A. This argument also holds for the lower edge of the 
conformal window where the magnetic theory is weak but for a larger number of flavors the 
size of k is unknown. 

When we add a chemical potential the dual quarks will get a negative mass-squared 
proportional to their baryon number, namely N^/Nc, in addition to the mass from (M). 
So altogether,^ 
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2{Nf/N,~l) - l^j 



(58) 



which leads to a critical chemical potential 

2 



k 



A 



(59) 



Thus for fi > the dual squarks will have tachyonic masses, signaling an instability. 
Because all the squarks get the same mass, we expect them to all get the same VEV x 
which, in order to satisfy the SU{N(.) D-term equations, is of the following form in terms 
of Nc X Nf matrices: 



V 



and 



X 



V 



(60) 



Such a VEV completely breaks the SU{Nc) gauge group, which means there is no residual 
U{1) that could be combined with U{1)b to give an unbroken U{1), hence we conclude that 
H exceeding the critical chemical potential will break baryon number. But we still need to 
check that such a value of fi is permitted within the validity of our approximations. 

We require Hc < A so that the strong-coupling results can be applied, which leads to an 
upper bound on m/A, 

NrA 




(61) 



*This formula assumes that the Kahler potential ^ d'^OZqq* q for the dual quarks is canonical {Zq = 1). 
If not, the normalization Zq appears in the formula by replacing k by Zqk. In fact, it is known that scales 
k, A, A are not physical, but rather depend on the normalization in Kahler potential This formula as 
well as the intuition fc ~ A ~ A make sense in the convention that the Kahler potential is canonical both 
for the electric quark in the UV limit and the magnetic quark in the IR limit. 
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The requirement that /Xc < m for stabihty in the UV leads to another constraint: 



The magnetic theory is free in the infrared, so we are justified in using the description in 
terms of dual quarks. Furthermore, the upper bound on m/A has a larger exponent than the 
lower bound this time, which means we need k/A > Nf/Nc — 1, which is a quite reasonable 
constraint since the right hand side is always less than 1/2. As argued above, we expect 
k and A to be comparable. This means that it is reasonable to expect k/A > Nf/Nc — 1; 
yielding a viable window where fi > fic exists and all our approximations remain valid. 
Therefore there is a critical chemical potential, of the form shown in the right diagram of 
Figure n above which baryon number is broken. 



7.1.2 The Conformal Window 

In this case we obtain the same formula for the dual quark mass as in Eq. (|58|) . The validity 
requirement fic < A also leads to the same constraint Eq. (jblf) . On the other hand, the 
requirement that fj,c < m for stability in the UV leads to three separate cases depending 
on the size of A'^^: 
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for Nf < 2Nc (63) 

for Nf > 2Nc (64) 
/ = 2iVe (65) 



We see that there are three situations. For Nf = 2Nc, f^c is directly proportional to m so we 
need A/k < 1. For more flavors, 2A''c < Nf < 3Nc, there are two different upper bounds, so 
it is easy to choose a small m that satisfies both. However, this region is the upper part of 
the conformal window where the magnetic theory is more strongly coupled, so our analysis 
in terms of the magnetic quarks is unlikely to be valid. Instead it is more appropriate to 
consider the electric theory where U{1)b is not dynamically broken by construction, since 
we require /i < m. 

The more interesting situation is when Nf < 2Nc- In this case m/A is bounded both 
from above and from below. For Nf = ^Nc the two bounds coincide, yielding jj.c = A, which 
doesn't leave any room for fi > within our approximations. For ^N^ < Nf < 2Nc the 
exponent of the lower bound is larger than the exponent of the upper bound, so for there 
to be any space between them, we need [N^/ Nc){k/ A) < 1, which means k/A < Nf/Nc — 1; 
which is always less than 1 for this range of flavors. While not impossible to satisfy, the 
naive expectation for Nf close to ^Nc is that A;/A ~ 1. 

We conclude that U{1)b remains unbroken in the conformal window. 
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7.2 SUSY-Breaking Mass 

We now wish to add soft masses for the electric squarks in the UV. Again we discuss the 
free magnetic phase and the conformal window separately. 

7.2.1 Free Magnetic Phase 

Since the soft mass modifies the electric theory we must modify the magnetic theory as 
well in order for the two theories to describe the same moduli space. In ^Uj it was shown 
by arguments similar to those of Appendix^ that the duality is maintained if soft masses 
are added to the magnetic quarks and the meson. In the deep IR these soft masses take 
the form 



mil 



3Nr - 2N 



f ~ 2 



and 



3Nr - 2N 



m 



<?.<? 



f ~ 2 




(66) 



where we have included the contribution due to the chemical potential. Note that the 
meson soft masses are positive, but as in the previous case, the dual quarks q and q get 
tachyonic masses. Again, since the magnetic theory is free in the IR, its quarks and meson 
are the appropriate degrees of freedom to analyze near the origin of moduli space. 

To analyze the symmetry breaking pattern first note that since the F-teim potential 
derived from Equation H56|l does not relate the expectation value of the dual quarks to the 
meson (or its cofactor), the meson VEV will always vanish. Representing the squark and 
anti-squark as two Nc x Nf matrices, the D-term equations can be satisfied with a VEV of 
the form 



Xl 



\ 



X2 



/ 



and 



Xl 



\ 



X2 



X'i 



(67) 



with \xi 



\Xi\ 



r where r is a common constant for all i. 



In order to find the vacuum we must minimize the rest of the potential (the non D- 
terms) with respect to Xi and the constant r. We can assume without loss of generality 
that r > 0. The potential including the soft terms is 



V 



\c2?'\Xi 



+r) + ^ml {2\xi\'^ +r 



(68) 



where is negative as in Equation (|66|) . This is minimized for Xj = with r going to 
infinity due to the negative mass squared. Since we have stabilized the theory for VEVs 
far away from the origin we expect r to stabilize at some finite value of order A. The 
instability of the origin exists even in the absence of a chemical potential since the dual 
squark masses in Equation (|66|) are negative. The case of Nf = ^Nc is an exception since 
then the masses of Equation H66|) vanish when // is set to zero, and r is a modulus. In that 
case an arbitrarily small chemical potential will destabilize the origin and break U{1)b- 
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Again we can only exclude the origin as a possible vacuum and point in the direction 
the theory will roll away from the origin. If the vacuum indeed lies along this direction 
the symmetry breaking pattern resembles the second flat direction for Nf = Nc + 'i with 
soft SUSY-breaking masses. Referring back to the quark VEVs in Equation (|67|) we see 
that this vacuum preserves the symmetry SU{Nc)r+c x SU{Nf — Nc)r x SU{Nf)L = 
SU{Nc)r+c X SU{Nc)r X SU{Nf)i^ ^. Negative r would correspond to the situation with 
left and right interchanged. We see here that the free magnetic phase exhibits some color- 
flavor locking in the sense that the (dual) color gauge group and a subgroup of the flavor 
symmetry are broken to the diagonal. 

7.2.2 The Conformal Windov^^ 

In it was shown that, unlike the supersymmetric masses that get enhanced, the soft 
masses in both the electric and the dual magnetic theories flow to zero in the IR. The rate at 
which this occurs is unknown and therefore we cannot determine whether the theory reaches 
the fixed point. Let us illustrate this qualitatively by focusing on the electric description. 
When the theory is weakly coupled (near the Banks-Zaks fixed point p^) the soft mass 
flows toward zero very slowly so we can approximate it to be constant at its bare value mo- 
The conformal symmetry will be broken at the scale of the soft mass. When the theory 
is more strongly coupled the soft mass fh{uj) is a more rapidly decreasing function of the 
scale uj. In this case the conformal symmetry will be broken at a scale that is equal to the 
soft mass at that scale to ~ fh{uj). In this case the IR theory will not be conformal nor 
supersymmetric, though the SUSY- breaking operator has been suppressed. 

Finally we may imagine that at yet stronger coupling the soft mass m(a;) flows down 
with scale faster than the scale uj itself. In this case the condition uj ~ m(a;) is never satisfied 
and the theory reaches a truly conformal fixed point. Note that even though supersymmetry 
was broken in the UV, it is restored in the IR. (See for an application of such a scenario 
as a solution for the hierarchy problem.) However, as was pointed out in |12j . since we 
cannot calculate the rate of approach to the fixed point we cannot determine whether this 
occurs. 

^From the qualitative picture just described we can get a rough idea about the size of 
the critical chemical potential for various numbers of flavors within the conformal window. 
Near Nf = ^Nc when the electric theory is strongest the IR soft mass is suppressed the 
most and a smaller chemical potential is needed to break U{1)b- Since the suppressed mass 
is obviously smaller than the UV mass, our assumption fi < ffiQ is satisfied and the U{\) 
breaking is indeed dynamical. If the last scenario of the previous paragraph occurs, baryon 
number is broken for an arbitrarily small /i. 

As we add flavors and move up in the conformal window the suppression weakens and 
the critical chemical potential rises. The naive expectation is that it will reach mo near 
Nf = 3Nc- Note that this matches the values of the critical chemical potential obtained 
outside of the conformal window: Hc = for Nf = ^Nc in Section 17.2.11 and fic = foi^ 
Nf > 3Nc where the theory is IR free, by the arguments of Section ^ The rising trend 

^Note that the SU{Nc) symmetry here is a flavor symmetry. 
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for the critical chemical potential is also consistent with the results for supersymmetric 
stabilizing masses in Section f7. II 

8 Conclusions 

Using known exact results for supersymmetric SU {N) gauge theories and soft masses, we 
have studied the effect of a baryon chemical potential ;U < A on the global symmetries of 
the ground state, in particular U{1)b- Our results are summarized in Table [3 



Nf 


Unbroken Global Symmetries 


SUSY Mass 


SUSY-Breaking Mass 


Nf<Nc 


SU{Nf)v X 


SUiNf)v X U{1)b 


Nf = N, 


SUiNf)v 


SUiNf)v 


Nf = Nc + l 


SU{Nf)v 


SU{Nf - 1) X SU{Nf) 
or SU{Nf - l)v 


Nc + 2<Nf < fiVc 


SU{Nf)v 


SU{N,)r+, X SU{N,)r X SU{Nf)L 
or {L ^ R) 


fiVe <Nf< 2N, 


SUiNf)v X U{1)b 


U{1)b broken 


2Nc <Nf < 3Nc 


SU{Nf)v X U{1)b 


U{1)b broken? 



Table 2: Unbroken subgroup of the original SU {Nf)L x SU{Nf)ji x U{1)b global symmetry 
for SQCD with the addition of either supersymmetric mass terms or soft SUSY-breaking 
masses. See the text for assumptions made in each case. 

We have found that for Nf < Nc a modified U{1)b remains unbroken at low energy, 
whereas for Nc < Nf < '^Nc there is a critical chemical potential above which baryon num- 
ber symmetry is broken due to a combination of strong gauge dynamics and the chemical 
potential. Notice that there is reasonable agreement between the SUSY-preserving and 
SUSY-breaking masses. Comparing with Tabled we find that U{1)b suffers the same fate 
as in non-supersymmetric QCD (Nc = 3) where it is unbroken in the 2SC phase {Nf = 2) 
but broken in the CFL phase {Nf = 3). For SQCD with Nf > 3 we also find interesting 
phases that include possible U{1)b breaking. For Nf = 4 we again find U{1)b breaking, in 
agreement with the non-supersymmetric case. Our results lend support to the results ob- 
tained using the variational approach, even though the latter are derived in a very different 
limit, namely S> A, and using very different methods. 
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A Harmonic Oscillator Example 

The two-dimensional isotropic harmonic oscillator is an explicit example which can be used 
to verify the formalism for treating the chemical potential. We begin with the Hamiltonian 
for the 2D harmonic oscillator and add a chemical potential for the "number" operator 

^ ~ ^ 2^^i^^ + y2) _ _ xpy) (69) 

which corresponds to the Lagrangian 

1 



L = —m 
2 



x^ + f- (u;2 - fi^){x^ + y^) - 2il{xy - yx)\ (70) 
where jl = ji/U. Writing z = ^{x + iy) puts this in the form 



L = m 



z* z - {u)^ - il^)z* z + iii{z* z - zz*) (71) 



which (apart from the factor of m) is the Lagrangian for a (0 + l)-dimensional charged 
scalar field with "mass" ui and nonzero chemical potential. 

Apart from the similarity between this system and the field theory that we are interested 
in, the virtue of this example is that we can compare three separate computations for the 
partition function. First, one can perform the sum over states in the grand canonical 
ensemble directly. This is most easily done by writing the Hamiltonian in terms of two 
sets of creation and annihilation operators, H = [Tiuj — ii)c^c +\ + {Tilo + iJ,)d}d + ^, where 
c = -^{dx + idy) and d = -^{o-x — icLy)- The partition function then becomes simply the 
product of two independent harmonic oscillators with frequencies offset by ib/x/^: 

2sinhte)2sinh^(^- ^ ^ 

The same result is obtained by computing the grand canonical partition function by an 
"imaginary time" path integral: 

Z = Tr e-'3(^-''^) = j VxVy e"^ ^^'^ dr[x-'+y''+{uj^-jX''){,x''W)+'iiKxy-yx)\ ^73^ 

where here x = drX. 
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We can also perform the real-time path integral computation of the quantum mechanical 
propagator, namely 



{Zf, tf\Zi, ti) 



Vz*Vz e 

iSc/n 



■7rihsin{uJ {tf — ti)) 



(74) 
(75) 



where 



Sr 



mu; 
sinwT 



l[ZfZi 



ZfZ*) sin fiT — {z%Zi + zjz*) cos /iT + {\zf\^ + kiP) cos ujT (76) 



is the action over the classical path, and T = tf — ti. Taking T = —ifif3 and Zf = Zi = z 
followed by integration over z gives the grand canonical partition function: 



Z 



dz{z\e 



1 



2sinh^^^P^2sinh^i^^^ 



(77) 



This matches the direct computation, verifying the form of the effective Lagrangian in 
Minkowski space. This exercise demonstrates how the /j,^ term in the Lagrangian is neces- 
sary to produce the correct result from a path integral calculation. 



B Exact Results for Soft Masses 

The following is a review of JJQ^ and [12j . We begin again with the SQCD Lagrangian 
renormalized at a UV scale ^c/y, 

Cuv = Jd^^T. Mf^uv)QlQr + J dPeS{iiuv) (W„W" + h.c.) , (78) 

where the index r labels the quark representations. SUSY breaking effects can be incor- 
porated by promoting the couplings Z and S to real and chiral superfields respectively as 
follows^*^ 



1 - e^Q^m^ 



Equation ()79() incorporates a universal soft mass m as well as a possible gluino mass m^. 
The theory possesses a spurious f/(l)A symmetry under which 

e-(^+^^) (80) 

and the gauge coupling transforms anomalously 

^°If there is no superpotential we can choose to absorb the and terms of Z into the gluino mass 
through the anomalous symmetry of Equations H8()|l and H81|l. 
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The main point of jlUj and |12j is that this symmetry is respected by the strong dynam- 
ics that lead to confinement. Since Z transforms as a gauge field under this C/(1)a, it's 
couplings to the low energy degrees of freedom are determined by this symmetry. This 
in turn determines the IR soft mass, which is the highest component of the superfield Z. 
Furthermore, any physical quantity, related to either the UV or the IR degrees of freedom, 
should be invariant under this symmetry and also RG invariant. It is therefore worthwhile 
to construct a \J(\)a and RG invariant object from S and Z since it will have physical 
meaning in both descriptions. The only such object is 



2JV 



K{z "0 Ah (82) 

where Ah = ^uv^~^'^^^^'^'^^'^^^'^ is the holomorphic dynamical scale and 60 = '^^c — is 
the one loop beta-function coefficient. 

Indeed, in Reference ^01 it was shown that for UV-free theories the various components 
of / are related simply to the physical dynamical scale and the bare soft terms 

[I]e=e=o = A' 

n n r IGvr^ rriA 
ln/g2 = hm — 2" 

[ln/]g2g-2 = lim m2 ^§3) 

bo fJ-uv-^oo 

To find the effects of soft masses in the low energy theory, all we need to do is write down 
the most general Kahler potential allowed by the symmetries, including the anomalous 
U{1)a, and determine the coupling of Z to the kinetic terms. Since we only need squark 
soft masses to stabilize the theory in the UV we will not add a gaugino mass. 

For example, consider the case with Nf = Nc+1- The Kahler potential can be expanded 
around the origin and is constrained by the symmetries and RG invariance to be 

K = CM—TvM^M + CBjj^Bje^^^^Bi + c^-^Bje-^^^^ Bi + ■■■ (84) 

where Vb is the background U{1)b gauge field. Note that since this is an expansion around 
the origin, the potential we will derive will only be valid for VEVs close to the origin that 
are much smaller than A. 

The soft masses for the mesons and baryons at some IR scale fijR will in general depend 
on the SUSY breaking VEVs of the spurion fields Z, I and Vb as well as 0{1) numerical 
coefficients cm, cb and c^. However, in ^] it was shown that due to the IR-freedom of 
the theory the soft masses don't depend on the coefficients cm-, cb, or in the deep IR, 
IJ-IR — *■ 0. Using the components of Z in Equation ((7^ . / in Equation (|55|) and Vb in 
Equation and performing the d'^O integral gives the soft masses for the canonically 
normalized mesons and baryons. Altogether we have: 

r^M = _ ^^'^ '^B,B = 2N - ■ (85) 

Notice that for Nc'>2> the meson mass-squared is positive, but the baryons are tachyonic. 
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